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SUMMARY 


The aerodynamic energy method is used to synthesize control laws for NASA’s drone for aerodynamic and 
structural testing-aerodynamic research wing 1 (DAST-ARW1) mathematical model. The performance of these 
control laws in terms of closed-loop flutter dynamic pressure, control surface activity, and robustness is compared 
with other control laws that appear in the literature and that relate to the same model. A control law synthesis 
technique that makes use of the return difference singular values is developed in the present work. It is based on the 
aerodynamic energy approach and is shown to yield results that are superior to those results given in the literature 
and are based on optimal control theory. Nyquist plots are presented, together with a short discussion regarding the 
relative merits of the minimum singular value as a measure of robustness as compared with the more traditional 
measure involving phase and gain margins. 

INTRODUCTION 

In recent years, a considerable amount of work has been done in the field of flutter suppression using active con- 
trols (Gangsaas and Uy-Loi, 1979; Mahesh et al., 1979; Mahesh et al., 1981; Mukhopadhyay, 1987; Mukhopadhyay 
et al., 1981, 1982; Mukhopadhyay and Newsom, 1982; Newsom, 1979; Newsom et al., 1983). The common denom- 
inator for most of these works is that they are based on optimal control theory. This requires the formulation of the 
aeroservoelastic equations in state-space form that results in equations involving many states. A realistic design, for 
example, with 10 structural modes, may lead to equations with approximately 60-70 states associated partly with 
the structural modes and actuator dynamics, and partly with the augmented states caused by the aerodynamic lag 
terms representation. 

Optimal control theory, based on the standard linear quadratic Gaussian (LQG) solution, yields an optimal feed- 
back control law of the same high order as the plant. Also, this optimal control law is invariably sensitive to modeling 
errors, has poor stability margins, and is often too complex to implement in a flight computer because of its high 
order. The robustness of the optimal control law based on the LQG design techniques can be considerably improved 
(Doyle and Stein, 1979) by introducing a fictitious input noise. However, the resulting LQG controller remains the 
same high order as the plant. Also, the minimization of the performance index involving the plant states cannot be 
readily related to the root mean square (rms) responses, since these former responses contain fictitious noise effects 
introduced into the system. As a result, a multitude of methods have been suggested for both reducing the order of 
the resulting control law and for maintaining the robustness improvement obtained by the introduction of fictitious 
noise. 

Methods suggested for obtaining reduced order control laws are derived using optimal control theory with trun- 
cation (Gangsaas and Uy-Loi, 1979; Mukhopadhyay et al., 1981), residualization (Mahesh et al., 1979), or transfer 
function matching (Newsom, 1979) methods. In the truncation method (Gangsaas and Uy-Loi, 1979), the low-order 
controller is obtained by partitioning the optimal full-state feedback and Kalman estimator gain matrices and re- 
taining only the optimal controller part corresponding to the key states and their first derivatives. Reoptimizing the 
preceding truncated control law using a numerical optimization algorithm (Mukhopadhyay et al., 1981), together 
with the improved robustness method using fictitious input noise, yields improved results (compared with the unop- 
timized truncation results using no fictitious noise). In the residualization method (Mahesh et al., 1979), a first-order 
correction is added to the truncation method (Gangsaas and Uy-Loi, 1979) by retaining the static part of the con- 
troller corresponding to the rest of the states. In the transfer-function matching method (Newsom, 1979), the optimal 
high-order control law is approximated by a low-order control law, the order of which is arbitrarily determined. The 
parameters of the assumed low-order control law are determined using an optimization algorithm that minimizes the 
deviations from the optimal control law over a range of frequencies. 

The main virtue of the optimal control theory is in providing tools for synthesizing control laws that stabilize 
the system using full-order controllers. These control laws clearly depend on the relative values of the different 



matrix elements assumed while formulating the performance index, the level of input noise (if any) introduced, the 
final order of the desired control law, and the method chosen for the control law order reduction. Also, the more 
potent of the previously mentioned design methods (Mukhopadhyay et al„ 1981; Newsom, 1979) require numerical 
optimization techniques to achieve the final control law. Therefore, the design of a control law involves a long 
process by which arbitrary parameters Oike those assumed in the performance index) are varied, and the resulting 
control law’s performance is tested and compared with the desired performance. 

The aerodynamic energy approach (Nissim, 1971, 1977, 1978; Nissim and Abel, 1978; Nissim and Lottati, 
1980; Nissim et al„ 1976 ) permits the stabilization of fluttering systems in a manner insensitive to flight conditions, 
mode shapes, and so forth, while using low-order control laws. Robust, low-order control laws form the basis of 
the aerodynamic energy synthesis technique. Its potency could be readily improved and its performance evaluated 
and compared with other design techniques based on optimum control theory. This latter comparison is the prime 
objective of the report. 

The report contains a brief review of the aerodynamic energy concept, together with the most important results 
obtained to date. Different control law design techniques, based on the aerodynamic energy concept, are evaluated. 
The results are compared with those based on optimal control theory. Finally, an evaluation of the aerodynamic 
energy control law design technique is presented. 

NOMENCLATURE 

Abbreviations 


DAST-ARW1 

DTTF 

LDTTF 

L.E. 

LQG 

MSV 

rms 

SISO 

T.E. 


drone for aerodynamic and structural testing-aerodynamic research wing 1 

damping-type transfer function 

localized damping-type transfer function 

leading edge 

linear quadratic Gaussian 

minimum singular value 

root mean square 

single input-single output 

trailing edge 


Letter and Mathematical Symbols 


ctL,aT positive gains 

6 semichord length at the sensor’s section 

c reference chord 

G, H matrices in open-loop transfer function 

g gravitational acceleration 

h plunge displacement, positive down 

h a plunge displacement of sensor 

M Mach number 

Q dynamic pressure 
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Q d design dynamic pressure 

Q f flutter dynamic pressure 

R u fictitious input-noise intensity matrix 

Sf transfer function of structural filter 

s the Laplace transform 

W work done on surrounding per cycle of oscillation 

W s weighting scalar defined in equation (13) 

a torsional angle, positive nose up 

(3 leading-edge deflection, positive nose down 

7 variable defined in equation (15) 

8 trailing-edge deflection, positive tail down 

8 C rotational command for trailing-edge deflection 

C damping coefficient 

\ aerodynamic eigenvalue 

cr singular value of return difference 

a minimum singular value of return difference 

desired minimum singular value 
u) frequency of oscillation 

uj n natural frequency of second-order pole 

Subscripts 

F flutter 

inb inboard driving sensor 

L leading-edge control 

max largest 

min smallest 

outb outboard driving sensor 

T trailing-edge control 

REVIEW OF AERODYNAMIC ENERGY CONCEPT 

General Approach 

The aerodynamic energy concept uses the work W done by the system on its surrounding per cycle of oscillation 
as a basis for determining the system’s stability. If the system does positive work per cycle of oscillation ( W > 0), 
the system is dissipating energy and is stable. But, if the system does negative work per cycle of oscillation ( W < 0 ), 
the system is absorbing energy and is unstable. At flutter, no exchange of energy takes place because the system is 
neutrally stable, and therefore W = 0 . 
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The mathematical expression developed (Nissim, 1971, 1977) to evaluate W shows that the aircraft s inertial 
and elastic forces have no contributions to W, except for control surface mass unbalance inertial forces where 
control surfaces are left unbalanced. When control surfaces are left unbalanced, the aerodynamic forces are the 
sole contributors to the energy exchange per cycle of oscillation. This result is not surprising because it is well 
known that dynamic systems involving only symmetrical mass and stiffness matrices always oscillate in neutral 
equilibrium and can never diverge in an oscillatory manner. Therefore, any flutter suppression effects from active 
controls are obtained by way of aerodynamic matrix changes made to the aerodynamic matrices using control surface 
aerodynamic contributions. In Nissim (1971, 1977), the expression for W is further reduced to a principal quadratic 
form involving quadratic coefficients and the square of the responses only (without cross products of responses). If 
these quadratic coefficients (denoted by Vs in Nissim (1971, 1977) and defined as aerodynamic eigenvalues) can be 
modified by an appropriate control law so that all become positive, then the system must be stable and dissipative. 
The square of the responses are always positive, thus yielding positive values for W regardless of the specific value 
of the responses. 

In an attempt to obtain general results and avoid treating specific systems, the Vs representing the previously 
mentioned coefficients for the principal quadratic expression for W were computed numerically (Nissim, 1971, 
1977) for a two-dimensional wing strip. Their values were modified through the use of generalized control laws 
to render all the Vs positive over a wide range of reduced frequencies and Mach numbers. The two-dimensional 
strip was allowed to plunge and pitch, and two control surfaces were permitted to be activated — a leading-edge 
(L.E.) and a trailing-edge (T.E.) control surface — each spanning 20 percent of the chord. The principal quadratic 
expression for W involves only two coefficients; that is, two X s. The smallest of these two X s is referred to as X TO j n 
and the largest is referred to as X raoI . The results obtained for the two-dimensional wing can be applied to the three- 
dimensional wings based on the following assumptions; The basic motion of each strip in the three-dimensional 
wing must involve plunge and pitch motions, and the ratios between the different aerodynamic forces acting on 
the three-dimensional strip must be the same as the ratios for the two-dimensional strip (much the same way as in 
Glauert’s lifting line theory). The main results pertaining to the two-dimensional strip (Nissim, 1971, 1977) are 
presented in the following section. 

Review of Aerodynamic Energy Two-Dimensional Strip Results 


The main results produced by the aerodynamic energy concept (Nissim, 1971, 1977) are summarized by the 
following points; 

1. An L.E. alone has very little effect on either X min or X rooi . 

2. A T.E. alone has little effect on X min and a large effect on X raai . With a proper control law, X moI can be made 
very large and positive, while X min is generally negative but can be brought to approach zero. 

3. An L.E.-T.E. system greatly affects both X TO ,-„ and X moi , with X TOin affected largely by the L.E. and X moI 
affected largely by the T.E. Both X m , n and X m0I can be made positive and assume large values. 

4. The most effective location on the strip for a single sensor activating the T.E. control is around the 
65-70-percent chord location, with a slight variation with subsonic Mach number used. 

5. The most effective location within the strip for a single sensor activating the L.E. control is around the 
-20-percent chord location (20-percent chord ahead of the strip’s L.E.). Since this location is impractical 
for a single sensor, two sensors need to be used to project the motion to the — 20-pcrcent chord location. 
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6. Denoting the motion of the T.E. and L.E. sensors located as described in items 4 and 5 by h a< x and h a< i, 
respectively, the most effective control law for all control systems (that is, T.E. alone and combined L.E.-T.E. 
systems) is where the T.E. control deflection is in phase with the velocity of its driving sensor (h a> x), and the 
L.E. control deflection is in antiphase with the velocity of its driving sensor (— 

Review of Specific Control Laws Obtained Using the Aerodynamic Energy Concept 

Using the notation presented in the following sketch, with h denoting the plunge displacement at the 30-percent 
chord location and with the arrows denoting the positive displacement and rotations, the fundamental aerodynamic 


Undisturbed position 



energy control law derived in Nissim (1977), and denoted as the damping-type transfer function (DTTF), is given 
by the following expression 



h/b 1 + T a L 
a J ax 


-4 4 ] J A/6 

4 3.2 J { d 


( 1 ) 


where ai and or are positive gains (or zero) determined while designing a specific system. For a T.E. alone system, 
ctL = 0 and only the equation associated with 8 need be considered. Also, when the gains ai and or are sufficiently 
large, the matrix associated with the displacements can be discarded. These coefficients affect the Vs very slightly, 
and use of the following simplified control law can be made: 


0 

8 


QL 

ax 


-4 4 

4 3.2 


(*) 


( 2 ) 


Equation (2) forms the basis for the results presented in the previous section under points 4, 5, and 6. Note that 
by placing the sensors at the locations described in points 4 and 5, the ratios between the two numerical columns in 
equation (2) are maintained. These ratios are important since multiplying factors can be lumped with the constants ai 
and ar (Nissim, 1977). The control law presented in equation (2) is of a basic nature and is shown to be approximated 
by the different robust control laws derived by a variety of design methods. However, this control law is impractical 
in its present form since it does not permit the introduction of damping over a limited range of frequencies. Instead, 
damping is introduced over an uncontrolled range of frequencies, with possible penalties paid in terms of unnecessary 
control surface activity. 

Also, equation (2) is impractical when considering the effects of actuator dynamics. The energy approach does 
not rule out the use of aerodynamic stiffness terms to change the response of the system for increase of responses 
associated with A TOOX . Thus the system is made more dissipative. This is true for a T.E. alone control system where 
\ m i n is small and not necessarily positive. Aerodynamic stiffness terms can be introduced through the terms propor- 
tional to the displacements, as in equation (1). These arguments led to the formulation of the localized damping-type 
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transfer function (LDTTF) (Nissim, 1977) summarized by the following control law that is a counterpart to the con- 
trol law presented in equation (2): 


0 

6 




s z + 2w ni iCz--s + 


a Tt 


s l + 


-4 4 

4 2.8 



(3) 


Again, a i and ax are positive gains, and s is the Laplace variable. The second-order transfer functions in the diagonal 
matrix in equation (3) yield terms proportional to s at resonance (that is, s = JWn,L or s = iuyr). Or, in other words, 
these functions yield control surface deflections in phase and antiphase with the rates of sensor displacements at 
the locations specified in points 4 and 5, much in accordance with the DTTF presented by equation (2). When at 
frequencies smaller or larger than the resonance frequencies of the second-order poles, these transfer functions also 
have real components (in quadrature with the velocity), and lead to aerodynamic stiffness terms in addition to the 
aerodynamic damping terms. The maximum amount of aerodynamic damping is introduced around the resonance 
frequencies of the second-order poles. These frequencies must be assigned numerical values around the flutter 
frequency. However, before summarizing the design procedure based on the foregoing results, it is stressed that the 
LDTTF is only one of the many transfer functions that can produce control deflections with large components in 
phase with h 3tT and in antiphase with h StL over a wide frequency range spanning the flutter frequency. Thus the 
fundamental results are presented by equation (2); that is, by the DTTF. Practical implementations can assume a 
very wide variety of forms , only one of which is given by the LDTTF. Additional examples of transfer functions 
yielding large components of control deflections in quadrature with the amplitudes of oscillation, as warranted by 
equation (2), are obtained by multiplying the LDTTF by lead-lag terms, or by considering the washout-type transfer 
function of the form s/(s+ a). In the following section, the design procedure developed in Nissim and Abel (1978), 
based on the LDTTF, is briefly presented. It is followed by the synthesis of control laws based on the foregoing 
results and on new techniques using a specific flutter example that includes actuator dynamics. 


Review of Control Law Synthesis Technique 

The idea behind the design technique developed in Nissim and Abel (1978) is to come up with a design that not 
only meets the stability requirements, but also yields minimum control surface activity. This is achieved by mini- 
mizing the response of the control surfaces to a gust input. In this manner the gust becomes a driver in establishing 
the closed-loop parameters, while maintaining stability and minimizing control surface RMS activity. The design 
procedure involves the following steps: 

1. Choose initial values for the free parameters ( o’s, w n ’s, and C’s ) in the LDTTF that stabilize the system. 

2. Minimize the control activity by optimizing the free parameters at the highest desired dynamic pressure and 
the highest subsonic Mach number, while constraining lo^l and w n> r to lie between 0.65 u;? and 1.45 uip 
where uj p denotes the flutter frequency. 

3. Determine control surface activity over a range of flight conditions (lower Mach number and lower dynamic 
pressures, and possibly other flight configurations). If considerable increase in control surface activity is 
observed, repeat step 2 with a narrowed range for uinj, and lj^x- 

The constraints on the values of and uj^x are intended to ensure that stabilization is achieved by changing 
the Vs, thus assuring robustness irrespective of the responses. Changing the responses of the system with little effect 
on the Vs may make it sensitive to changes in responses brought about by small changes in flight conditions or by 
small modeling errors. 
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DAST-ARW1 MODEL 


The objective of this work is to evaluate the aerodynamic energy design method by comparing the performance 
of its resulting control laws against those obtained by other methods based on the different forms of the optimal 
control theory. The mathematical model for NASA’s drone for aerodynamic and structural testing-aerodynamic 
research wing 1 (DAST-ARW1) was chosen for this purpose because it served as a basis for the development of the 
different design techniques based on the optimal control theory. However, the complete mathematical model was no 
longer available. It had to be reconstructed using generalized mass and stiffness matrices, together with the mode 
shapes associated with the model’s different degrees of freedom. The aerodynamic matrices were computed using 
the available mode shapes supplied by NASA. Therefore, the first step in the present work is to validate the model 
against the published results for the open- and closed-loop systems, bearing in mind that some variations might exist 
as a result of possible differences in the aircraft’s aerodynamic modeling. 

Description 


The DAST-ARW1 model (Newsom ct al., 1980) is a dynamically sealed representation of a transport-type re- 
search wing and was sealed to flutter within the operational limits of the Langley Transonic Dynamics Tunnel. The 
model was equipped with a hydraulically actuated T.E. control surface located between the 76.3 percent and 89.3 per- 
cent semispan stations. Its chord measured 20 percent of the local wing chord. The model’s geometry is given in 
figure 1. 

The mathematical model used for the acroclastic analysis and acroscrvoclastic synthesis comprises 10 clastic 
symmetrical modes. Its natural frequencies span the range of 5.23 to 118.15 Hz. 

The unsteady aerodynamic forces for the wing and control surface were computed at Mach 0.9 for different 
values of reduced frequencies using doublet lattice aerodynamics. Unsteady aerodynamic forces were calculated for 
the 10 structural modes, the control surface rotation, and a sinusoidal gust. 

The T.E. control surface, driven by an clcctrohydraulic scrvoactuator, was capable of approximately ±14° rota- 
tion with rates of up to 820 dcg/scc. The actuator transfer function is given by 


6 _ 1 .915 x 10 7 

8~ c ~ (s+ 214)(s 2 + 179 ,4s + 8.945 x 10 4 ) dCg/ Cg 

where 5 C denotes the rotational command. 


(4) 


Validation 


Newsom ct al. (1980) form the basis for most of the comparisons between the present and the previously pub- 
lished performances. A more limited use is made of the results published in Mukhopadhyay ct al. (1981). Most of 
these latter results relate to reduced-order systems while attempting to synthesize relatively low-order control laws. 

Figures 2 and 3 compare the root locus plots of the open-loop system. The plots arc identical, except for a small 
difference in the flutter dynamic pressure (Qf), with Newsom et al. (1980) yielding Qp =105 lb/ft 2 compared with 
Qf - 99 lb/ft 2 obtained here. There is also a small difference between the flutter frequencies (u;^), with Newsom 
ct al. (1980) yielding ujp = 51 .6 rad/scc compared with ujp = 50 rad/scc obtained here. 

For the closed-loop system, comparisons arc made between the closed-loop values of Q p and u ;/?, Nyquist plots, 
stability margins, and control surface activities. 
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Figure 4 shows the closed-loop root locus plot obtained using the current model and the control law given in 
Newsom et al. (1980) by equation (9); that is 


Jc_ 

houtb 


2214 a 

s 2 + 30 .79 a + 14,692' 

a 2 + 47 .37 a + 72,436' 

deg 

(s + 10) (s + 1) 

a 2 + 572 .6 a + 88,578 

a 2 + 568 .6 a + 86,972 

9 


(5) 


where houtb refers to the accelerometer output located outboard of the T.E. control at the 60-percent chord location 
(fig. 1). In addition, a structural filter was added to the actuator in Newsom et al. (1980) to attenuate the high 
frequencies. This filter is denoted by Sf and is given by 


Sf = 


(628) 2 
(a + 628) 2 


( 6 ) 


Figure 5 shows the closed-loop root locus plot reproduced from Newsom et al. (1980). Note that figures 4 and 
5 are similar, with the present work yielding Qf = 174 lb/ft 2 and wf = 49 rad/sec compared with the values of 
Qf = 188 lb/ft 2 and wf = 44 rad/sec given by Newsom et al. (1980). The percentage differences between the 
values of Qf obtained in this report and those obtained in Newsom et al. (1980) are approximately the same for both 
the open-loop and closed-loop cases. 

The design point in Newsom et al. (1980) relates to a dynamic pressure of 44 percent above the open-loop flutter 
dynamic pressure. Since Qf = 105 lb/ft 2 in Newsom et al. (1980), the design dynamic pressure, denoted by Qd, 
was assigned the value of Qd = 151 lb/ft 2 . However, since Qf = 99 lb/ft 2 in the present work, the design dynamic 
pressure was 44 percent above this value, yielding Qd = 142 .6 lb/ft 2 . Figures 6 and 7 show a comparison between 
the Nyquist plots obtained in this report and in Newsom et al. (1980). These figures are very similar, with figure 6 
relating to the results obtained in the present work and yielding gain margins of —5 .34 dB and 14.3 dB, compared 
with —6 .0 dB and 16.5 dB obtained in Newsom et al. (1980). Similarly, the phase margins obtained are given by 
—48° and 47°, compared with the Newsom et al. (1980) values of —52° and 52°. The minimum singular value 
(MSV) of the return difference computed in this work is equal to MSV = 0.61. Because Newsom et al. (1980) does 
not present any values of the MSV’s, these values were extracted graphically from the Nyquist plots (explained in 
detail later in the report), yielding for this case the value of MSV = 0.64. All of the results, together with results 
pertaining to control surface activities and additional results presented in the following sections, are summarized in 
table 1. 

A comparison is also made between the results obtained in this report and those reported in Mukhopadhyay et al. 
(1981), using their control law, in equation (30): 

S c _ 1 939 . 4( a + 24 .74) (a 2 + 87 .63s + 13,806) deg 

(s+ 3.864)(s+ 3270)(s 2 + 20.97s+ 1423) g { * 

The same actuator used in Newsom et al. (1980) was also used in this case. However, no use was made of the 
structural filter given by equation (6). In this instance, the design point in Mukhopadhyay et al. (1981) was chosen 
to be at Qd = 160 lb/ft 2 . Considering the differences between Newsom et al. (1980), Mukhopadhyay et al. (1981), 
and the present work regarding the open-loop flutter dynamic pressures, the equivalent value of Q d for comparison 
purposes is Qd = 151 lb/ft 2 (= 160 x ^ ) . The results pertaining to the control law given by equation (7) (or 
eq. (30), Mukhopadhyay et al. (1981)), are summarized in table 1 and shown in figures 8-11. Note that there is good 
agreement between the results computed here and those reported in Mukhopadhyay et al. (1981). 

Since these results form the basis for comparisons with the control laws derived in this report, mention is made 
of the methods employed for synthesis of these control laws. The control laws derived in Newsom et al. (1980) are 
based on optimal control theory, followed by control law order reduction using transfer-function matching leading 
to a sixth-order control law. 
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The control law derived in Mukhopadhyay et al. (1981) used a reduced-order model for the DAST— ARW1, 
yielding a 25th-order plant instead of the original 65th-order plant (full order). Also, truncation and reoptimization 
were used to derive the control law, together with the addition of fictitious noise at the input to improve robustness. 
Finally, the reduced-order control law thus obtained was evaluated using the full-order model of the plant. 


DEVELOPMENT AND EVALUATION OF IMPROVED CONTROL LAW 
SYNTHESIS TECHNIQUE 

In the following, control laws are synthesized in three cases. In the first case, actuator dynamics are ignored, 
and a procedure is developed to replace the relatively loose optimization constraints imposed on uinj' by a more 
definite requirement, such as the increase of the minimum singular value over a wide range of frequencies. The 
manifestation of the MSV on the Nyquist plots and its relation to phase and gain margins is studied to define an 
improved-design procedure. In the second case, actuator dynamics are introduced, together with the structural filter, 
using the design procedure defined in the first case. The results obtained in the second case are compared with those 
given in Newsom et al. (1980) and presented in table 1 . In the third and final case, actuator dynamics are considered. 
However, the structural filter introduced in Newsom et al. (1980) is ignored, much the same as in Mukhopadhyay et 
al. (1981). The results obtained during the third case are compared with those given in Mukhopadhyay et al. (1981) 
and summarized in table 1. In all three cases, control laws are synthesized using the outboard accelerometer. The 
results are compared with the equivalent control laws synthesized using the inboard accelerometer that is warranted 
by the aerodynamic energy approach. Finally, the results obtained by the aerodynamic energy approach are evaluated 
and conclusions are drawn. 


Case 1 — Synthesis Technique and Presentation of Results for Systems With No Actuator Dynamics 


As stated earlier, the DAST-ARW1 mathematical model consists of a wing with only a T.E. control surface. 
Thus, equation (3) can be simplified and brought to the form 


6 = 


4ot'S 2 


a 2 + 2w n> rCr« + 


r°- 7 “) 


( 8 ) 


Because h relates to the 30-percent chord point, the expression (h + 0 .7 6a ) relates to a point at the 65-percent chord 
location. Therefore, equation (8) can be written as 


6 = 


4 (ITS 2 


fio.6Sc 


S 2 + 2w ni rCrs + b 


( 9 ) 


or 


6 = 


4ars 2 


hinb 


s 2 + 2uvr0rs+ w 2 r b 


( 10 ) 


Equation (10) is valid for the inboard accelerometer mentioned earlier. For the outboard accelerometer located on 
the 60-percent chord point, this location can be written by 


ho.6 c = fio,3c+ 0.66a 


( 11 ) 


Equation (11) leads to an equation similar to the one given by equation (3), but with the value of 2.8 replaced by 2.4 
(0.6 x 4). This is a suboptimal location because of the change in value and the sensor’s location just outside the 
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strip containing the T.E. control surface, instead of along its center span (fig. 1). For the outboard accelerometer, the 
following expression is used: 

^ _ 4 a T s 2 hp ,6 c 

s 2 + + 6 

( 12 ) 

or 

4 h ou tb 

0 = — r 5 — V 1 

+ 2w <r Cr$ + b 

Following the aerodynamic energy control law synthesis procedure, the values of ar » » and C t are evaluated 
both for the inboard and the outboard accelerometers. The target function used involves the rms value of the rate 
of the control surface response 8^ to a 1 -ft/sec gust input. A compounded target function can be defined that also 
includes the control rms response in addition to the control rms rate; however, this possibility was not used. The 
Von-Karman gust spectrum is used throughout this report, with a characteristic length of 100 ft, to simulate the 
distribution and intensity of turbulence. 

All the results obtained during case 1 are summarized in table 2 and the control laws are presented in appendixes A 
and B. Appendix A relates to control laws prefixed by the letter A and synthesized while using the inboard accelerom- 
eter as a driving sensor. Appendix B relates to equivalent control laws prefixed by the letter B and synthesized while 
using the outboard accelerometer as a driving sensor. 

The first optimization runs were intended to test the effects of the optimization constraints imposed onm^. To 
show the effects of the constraints on c , the value of was constrained to lie within the range 2 < < 

70 rad/sec. In the second example, was constrained to lie within a narrower range defined by 35 < < 

70 rad/sec with the aerodynamic energy synthesis technique. Note that only the lower bound was changed, since it 
is this bound that produces an active constraint. 

Control laws A1 and B1 relate to the case where is permitted to go as far down as = 2 , where it hits 
the constraint. With u - 2 rad/sec and with the flutter frequency ojp = 50 rad/sec, it is unrealistic to expect large 
components (approximately the flutter frequency) of control surface deflections that are in phase with the velocity of 
theT.E. driving sensor. Therefore, although the control surface rates yield values of approximately 174 deg/sec, 

the phase margins obtained are poor (approximately 20° to 24°). The MSV’s obtained assume the values of 0.31 for 
the outboard accelerometer and 0.35 for the inboard accelerometer. At this stage the design procedure based on the 
aerodynamic eneigy concept, as presented in Nissim and Abel (1978) and described earlier in this report (see point 2 
of the reviewed synthesis technique), does not permit such a wide deviation of from the flutter frequency. 

Control laws A2 and B2 relate to the case where the lower bound of is increased to 35 rad/sec in adher- 
ence with the design procedure presented in Nissim and Abel (1978). The resulting control laws (A2 and B2) have 
smaller gains ax and smaller values of as compared with those associated with control laws A1 and Bl. Thus, 
control laws A2 and B2 yield most of their responses in phase with the sensor’s velocity over a narrower range of 
frequencies, with the response peak located much closer to the flutter frequency. These control laws result in sharp 
increases in the MSV’s, accompanied by large increases in the positive phase margins. Therefore, the inboard ac- 
celerometer yields MSV = 0 .59 , with phase margins of —65 ° and 50° and gain margins of —6 .61 and +oo dB. The 
outboard accelerometer yields MSV = 0 .56 , with phase margins of —68° and 45° and gain margins of —7 .31 and 
10 .3 dB. This increase in robustness is consistent with the predictions made by the aerodynamic energy approach. 
It is achieved at the expense of a considerable increase in the minimum control surface rms rates that reach values 
of approximately 207 to 211 deg/sec. 
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To improve the synthesis technique, an attempt was made to incorporate the minimum singular value into the 
target function J so that 

J= 6™ + W s (s_-o D ) 2 (14) 

where £ denotes the MSV over the whole frequency range, and denotes the desired value of £. The factor 
W s represents the relative weight given to the MSV in the target function. The optimizer was modified to handle 
discontinuities. Control laws A3 and B3 were obtained using W s = 1000 , 35 < < 70 rad/sec, and a p = 0 .9 

for control law A3 and a_ D = 0 .8 for control law B3. The results obtained are shown in figures 12 and 13 and in 
table 2. The main effect of introducing the MSV into the target functions is to further reduce the gains ap- Also, 
although the MS V’s increase so that £ = 0.65 for control law A3 and £ = 0.64 for control law B3, the negative gain 
margins are reduced in magnitude compared with control laws A2 and B2. It is difficult to observe improvements in 
phase margins in control law A3 when compared with control law A2 even if the sign of the margins is considered 
unimportant. Yet there is a definite improvement in the value of a. This point is discussed in more detail later in the 
report. The values of the MSV’s and the values of the phase and gain margins produced by control laws A2, A3, 
B2, and B3 are very close to those presented in table 1, based on the optimum control theory, and even exceed them 
in some cases. 

The performance of control laws A3 and B3 at dynamic pressures Q below the design dynamic pressure Qd 
are tested at Q = 134 lb/ft 2 and Q - 124 lb/ft 2 . The results are summarized in table 2. The decrease in dynamic 
pressure is accompanied by a slight increase in the MSV’s and a significant increase in both phase and gain margins. 
Thus, discrepancies between the MSV’s and the phase and gain margins are again observed. As expected, control 
surface activity reduces significantly with the decrease in Q. 

The performance of control laws A3 and B3 in flutter suppression is shown in the root locus plots in figures 14 
and 15. The flutter dynamic pressure is increased from Qp = 99 lb/ft 2 to Qp = 175 lb/ft 2 for control law A3 and 
to Qp = 178 lb/ft 2 for control law B3. 

The previously mentioned small changes in the MSV’s as the dynamic pressure was lowered, using control laws 
A3 and B3, prompted the next examples in this report. This involved the reoptimization of control laws A3 and B3 
assuming Qd = 124 lb/ft 2 . The results are presented in table 2 as control laws A4 and B4. There is a noticeable 
increase in the MSV’s leading to £ = 0.83 for control law A4 and £ = 0.82 for control law B4. The phase margins 
are also improved. The smallest improvements are obtained in the negative gain margins. In all cases the control 
surface activities decrease with Qp>. 

Finally, control law B3 was reoptimized by removing the lower bound on much below ~ 35 rad/sec, 
while keeping the target function in the form of equation (13) with W s = 1000 . The results obtained were identical 
to those of control law B3, thus indicating that the constraint is effectively nonactive. 

The results presented in table 2 and in figures 12-15 can be summarized by stating that the bounds imposed on 
u) nt T by the aerodynamic energy approach can be substituted by a target function of the form given by equation (13) 
involving minimum singular values and the rms deflection rates. Control surface activity is always reduced as Q 
and Q d both decrease. For a given Qd , small changes in g_ are observed as Q is reduced. However, large changes in 
a are observed as Qd is reduced, thus indicating a limit to robustness obtained for any chosen Qp>. This robustness 
decreases as Qd increases. If both robustness and Qd need to be increased, additional or larger control surfaces must 
be considered. The term robustness is used loosely since improvements in phase and gain margins do not always 
correlate with improvements in MSV’s. However, the target function of the form of equation (13) is adopted for the 
remainder of this work as a substitute for the constraints imposed on 

In all cases, robustness is increased at the expense of increased control surface activities. This manifests itself 
in terms of control laws with smaller gains of ax and smaller values of 
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Case 2 — Control Laws Synthesized With Actuator and With Structural Filter 


Definition of the Form of the Required Control Law 

The general form of control laws treated thus far can be written as 

8 = 1 3L (15) 

(s 2 + & 

where h s denotes either the outboard or the inboard sensor. The introduction of an actuator and a structural filter 
clearly affects the control law so that the actual control deflections will no longer have a large inphase component, 
with the velocity of the driving sensor over a wide range of frequencies spanning the flutter frequency. Also, if 
sufficient care is not exercised, the control law— actuator combination may lead to control deflections with large 
components in antiphase, with the velocity of the driving sensor at relatively low frequencies leading to instability 
of otherwise stable modes. 

The basic approach adopted in this report assumes equation (14) to be the desired control deflection-sensor 
movement relation. The actuator is compensated for as much as possible, with more weight given to the lower 
frequency range. Following this initial compensation, the resulting control law is optimized in much the same way 
as in case 1 . 

Starting first with the structural filter given by equation (6), note that the filter’s poles are high enough so that their 
effects on the low frequencies adjacent to the flutter frequency can be neglected with no need for compensation. The 
effects of the actuator given by equation (4) can be partly compensated and partly minimized by writing equation (14) 
in the following form while replacing 6 by 6 C : 

6 = iHJll Aar tf + W +rih (16) 

c 7173(3 + 74) (a 2 + 2w„ i rCr3 + w^y) 6 

The final control deflection-sensor movement relation is given by the following expression using equations (4), (6), 
and (15): 

1.915 x 10 7 74 (s + 71) 4 qy( a 2 + 72 a + 73 ) 

^ - (s+ 214)(s 2 + 179 .4 s+ 8.945 x 10 4 ) X 7173 (3 + 74) (s 2 + 2w n> yCys + w^y) 

(628) 2 h, m) 

X (s+ 628) 2 b 

Initially, if 71 = 214,72 = 179.4,73 = 8.945 x 10 4 ,and7 4 is assigned a high value such as 74 = 630, then 
the numerators in equation (16) cancel the actuator poles. All that remains is the relation shown in equation (14), 
together with additional poles that affect only the high-frequency range. These poles could be shifted to even higher 
values, if necessary. Equation (16) is clearly dependent on the form of the actuator transfer function. The desigr 
procedure involving the previous control law requires the initial assignment of numerical values to all the constants 
in equation (17). The 7’s are assigned the values stated earlier, with 71 slightly reduced to compensate for the phase 
lags introduced around the flutter frequency by the structural filter and by (s + 74). The parameters a T , Wn,y . and 
are initially assigned the values associated with control laws having no actuator, as those obtained in case 1 of 
this report. Once the system is stabilized at Q p following the assignment of the previous initial values, the response 
of the system to gust is used in much the same way as in case 1 to yield optimum values for all 7’s, ay, w^y, and 
£y. More parameters are now available to the designer during the optimization stage, with more freedoms for the 
system to reach a possibly better optimum. 

This manner of introducing the actuator is more consistent with the aerodynamic energy approach than is the 
Newsom et al. (1980) approach. All phase compensations in Newsom et al. (1980) were made through a single lead 
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term that restored the phase values only around the flutter frequency. Clearly, at frequencies other than the flutter 
frequency, such a method yields a control-sensor output relation that is far different from the one warranted by the 
aerodynamic energy method. 


Presentation of Results 

Control law A5 uses the inboard accelerometer and control law B5 uses the outboard accelerometer. These con- 
trol laws were obtained following the design procedure previously described using the parameters associated with 
control laws A3 and B3 designed in case 1. The resulting optimized control law parameters are presented in ap- 
pendixes A and B. The results associated with these control laws are presented in figures 16-23 and are summarized 
in table 3. 

In table 3, the MSV’s are slightly larger than those obtained for control laws A3 and B3, with g_ = 0 .67 for both 
control laws. Similar smaller improvements are also noted in the phase margins (—55° and 59° for control law A5 
and — 61 ° and 60° for control law B5). More significant improvements in mns control deflection rates occur, yielding 
firms = 204.3 deg/sec for control law A5 and fi^s = 206.8 deg/scc for control law B5. The Nyquist plots associated 
with these two control laws are shown in figures 16 and 17 at the design dynamic pressure Qd~ 142 .6 lb/ft 2 . 

Table 3 summarizes the performance of these control laws at dynamic pressure lower than Q £>, and the associated 
Nyquist plots are presented in figures 1 8-2 1 . The closed-loop root locus plots are shown in figures 22 and 23. Values 
for the MSV’s remain unchanged as Q is reduced to Q = 134 lb/ft 2 and Q = 124 lb/ft 2 . However, both the phase 
and gain margins show significant improvement as Q is reduced. Some clarification of requirements is necessary, 
and a decision must be made as to what is more important — phase and gain margins or minimum singular values. 
This conflict is discussed in a later section of the report. 

Finally, these improvements in robustness of control laws A5 and B5, compared with A3 and B3, are accom- 
panied not only by a decrease in control surface rms activity, but also by an increase in closed-loop flutter dynamic 
pressure. In figures 22 and 23, the flutter dynamic pressures associated with control laws A5 and B5 are given by 
Qp - 193 lb/ft 2 and Qp = 189 lb/ft 2 , respectively. These values are higher than the values obtained for control 
laws A3 and B3 designed for the case where no actuator and no structural filter are present. Thus, the introduc- 
tion of the actuator and filter into the design procedure does not degrade the resulting performance, but noticeably 
improves it. 

Case 3 — Results for Systems With Actuator Only, Without Structural Filter 

The structural filter was introduced in Newsom et al. (1980) based on engineering judgment only. However, the 
work in Mukhopadhyay et al. (1981) used the same mathematical model of the DAST-ARW1, but the authors chose 
not to use a structural filter. Therefore, the work described in case 2 is repeated in case 3, but without the structural 
filter defined by equation (6). 

The results obtained are summarized in table 4 and are presented in part in figures 24 and 25, with control law 
parameters presented in appendixes A and B. Note that the removal of the structural filter results in both a reduction 
in control surface activity and in an increase in robustness, both in terms of MSV’s and in phase and gain margins. 
The MS V for control law A6 associated with the inboard accelerometer is cr = 0 .74 , and for control law B6 the MS V 
associated with the outboard accelerometer is g_ - 0 .79 . An improvement in phase and gain margins can also be 
noted, with control law A6 yielding phase margins of — 64° and 61°, and with gain margins of —5 .01 and 15.39 dB. 
Control law B6 yields phase margins of — 69 ° and 62 0 , with gain margins of —5 .3 and 14.33 dB. The control surface 
rms rate for both control laws is approximately 191 to 192 deg/sec. The performance of these two control laws at 
dynamic pressures below is summarized in table 4. As Q decreases, the MSV’s vary slightly, whereas both the 
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phase and gain margins improve considerably. Also, in agreement with all previous control laws, control surface 
activity decreases with Q. The flutter dynamic pressure is increased for both control laws to Qp = 206 lb/ft 2 . 

Finally, control law B6 was redesigned at a higher Q D using the outboard accelerometer only. This was done at 
Q d =151 lb/fl 2 to compare the results obtained in this report with those reported in Mukhopadhyay et al. (1981) 
at Qd = 160 lb/ft 2 . The ratio between these two values for Q d is identical to the ratio between the open-loop 
flutter dynamic pressures obtained in this report and in Mukhopadhyay et al. (1981). The results relate to control 
law B7, and are summarized in table 4 and presented in figures 24 and 25. Note that despite the previously mentioned 
increase in Q D , the MSV and the control surface deflection rate are only slightly reduced. 

Minimum Singular Values Compared With Phase and Gain Margins 

It is often stated that MSV’s lead to conservative estimates regarding bounds of phase and gain margins. There- 
fore, they should not serve as a means of actual design of robust control laws (too conservative). This investigation 
involves a single input-single output (SISO) system and is very convenient for the consideration of this issue because 
Nyquist plots are meaningful for SISO systems and their relation to the MSV’s can be readily considered. 

Consider the return difference (1 + GH) of a SISO system. The singular value a of this system at a given 
frequency is given by the absolute value of the scalar ( 1 + GH), whereas the Nyquist plot represents the values of 
the complex scalar GH at the different frequencies. Thus, one can write 

a = \(GH — ( — 1))| 


or 


a - R 


where R represents the radius vector joining any point from the Nyquist diagram to the —1 point, as shown in the 
following sketch. Clearly R is a function of the frequency w and, therefore, the minimum singular value a is equal to 
the minimum value of R connecting point —1 to the Nyquist curve, as with point C in the sketch. The MSV relates 
to a point where the combined phase and gain variations are closest to the unstable condition. 
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To demonstrate the problems which may arise with phase and gain margins, consider the following Nyquist 
diagram. Both the gain and phase margins are very good. However, a small variation in gain can drastically reduce 



Real, rad/sec 


9003 


Imaginary, 

rad/sec 


the phase matins. A Nyquist plot of this type is not atypical. An example similar to this sketch is given in figure 19. 
Other figures in this report can also support the argument. Similarly, small changes in phase margins can lead to 
large changes in gain margins, as illustrated in the following sketch. Figure 6 shows that this is possible, although in 
a less extreme form. This illustrates that there is no direct correlation between MSV’s and the phase and gain mar- 
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gins. On the basis of the results obtained in this report and the previous discussion, it can be stated that phase 
and gain margins do not present a proper criterion by which robustness can be measured. Future designs should 
relate to MSV’s only , and no attempts should be made during the design process to associate these values to ei- 
ther phase or gain margins. Therefore, in the remainder of the report robustness will be measured in terms of 
MSV’s only. 

Discussion of Results 

The results shown in table 2 relating to control laws designed with no actuators indicate that the constraints 
imposed on the values of and by the aerodynamic energy approach are necessary to ensure the robustness 
of the synthesized system. This point explains why the design based on the aerodynamic energy concept in Newsom 
et al. (1980) yielded relatively low robustness. Newsom ct al. (1980) imposed almost no constraints on the values 
of Once the constraints are imposed, yielding control laws A2 and B2, the MSV’s were fairly close to 
those obtained while trying to maximize a, yielding control laws A3 and B3. Also, in the range of values assumed 
for Qd> the values of a became significantly larger as Qd was decreased. However, for a given value of Q^, the 
value of a varied slightly as Q was reduced. This latter effect is shown in all of the tables. 

Some interesting results of this report relate to those obtained when actuator dynamics were introduced (see 
tables 3 and 4). The inclusion of actuator dynamics did not deteriorate the performance, but actually gave better 
results than those associated with no actuators. Table 3 shows that although the MSV’s are only slightly improved 
when compared with table 2, control surface activity is significantly reduced. Further improvements are also shown 
in table 4 with the structural filter removed. Again, the improvements encompass both the MSV’s and the control 
surface activities, with all minimum singular values assuming values above a = 0 .74 and with some MSV’s reaching 
as high a value as a = 0 .79 . The results presented in tables 3 and 4 exceeded those presented in table 1 that were 
obtained by various authors based on variations of optimal control theory. Higher order control laws introduce more 
freedoms through which robustness can be improved, together with reduced control surface activities. 

Examine whether the resulting control laws exhibit a large component of control deflection that is in phase with 
the velocity of the sensor. Plots of 8/h s as a function of frequency are shown in figures 26-33. For each control law, 
two plots are presented: one with the imaginary component of 8/h s that is in phase with the velocityat the sensor’s 
location when it assumes positive values, and one with the real component of 8/h s . Plots are presented for most 
of the control laws appearing in tables 1-4, with an additional control law taken from Mukhopadhyay et al. (1981), 
equation (29). The control law taken from this latter equation was derived using the truncation and reoptimization 
method, but without the robustness recovery obtained through the introduction of a fictitious noise at the input. 

The resulting control law yields poor phase and gain margins. Its 8fh s plots are shown in figure 28. In the future, 
control laws with low values of a have either smaller positive imaginary values of 8/h s or else the positive imaginary 
values span a more limited range of frequencies. This again is in accordance with the results of the aerodynamic 
energy approach. When 8/h s has a negative imaginary part (that is, control deflections are in antiphase with the 
velocity of the driving sensor), the T.E. active control system is likely to have a destabilizing effect. However, if this 
occurs at sufficiently high frequencies, the neighboring structural modes should be stiff enough and should remain 
stable. This is true provided the values associated with the negative imaginary parts of 8/h s are kept reasonably 
small. 

The simplified design procedure, based on the LDTTF with no actuators, yields results only slightly inferior to 
those obtained using actuators and compensating lead-lag terms. Thus, this simplified design procedure can be used 
during preliminary design stages, or in conjuction with large packages involving optimization of the structure along 
with its active control system. Final designs that include actuators can be left to the final stages after the structure 
assumes its final form. 


16 



Finally, this report shows that no distinct relative advantages can be attached to either sensor location. Nonethe- 
less, the inboard location is advantageous. The actuator reacts directly to the motion at the center of the strip where 
it is located, rather than to the motion of a remote point that bears no relation to the motion of the active strip. Addi- 
tional work involving other Mach numbers, or additional configurations, are needed to show the relative advantages 
of the inboard sensor location. 

The optimizer used in this report is a variation of Stewart s adaptation of the Davidon, Fletcher, Powell method 
(Fletcher and Powell, 1963; Stewart, 1967). The variation introduced permits the values of the optimization param- 
eters to be constrained without using penalty functions as described in Nissim and Abel (1978). 

Shortly after the introduction of equation (13) as a target function, the optimizer exhibited difficulties in con- 
verging. These difficulties were eventually traced to the discontinuous nature associated with the minimum singular 
values. For purposes of illustration, assume that (1) the term associated with the singular values is the dominant 
term in equation (13); (2) two singular values have the same minimum value (that is, a\ = oi = £, as shown in 



the sketch); and (3) a parameter change that increases any one of these singular values decreases the other one, and 
conversely. This introduces a discontinuity, making it impossible for the optimizer to make any progress. 

To counteract such a situation, the optimizer was modified so that it tested for additional singular values within a 
chosen small radius from the value of o_ at frequencies that arc not immediately adjacent to the frequency associated 
with q_. When such points were found, the derivatives with respect to the design parameters were evaluated by both 
forward and backward differencing. If the derivatives evaluated by forward differencing were found to change their 
sign when computed by backward differencing, then those derivatives were assigned the value of zero. Once these 
modifications were introduced, no convergence problems were encountered. 

Summary of Proposed Control Law Synthesis Technique 

The following control law synthesis procedure is suggested based on the foregoing discussion. However, varia- 
tions to this procedure can be made once sufficient experience is gained using this method. 

1. Choose the localized damping-type transfer function with no actuators. 

2. Choose initial values for q.t< ^n.T > and ^ ( or ai, ^n,L > and Cl) for ah active control surfaces to stabilize the 
system at Qd and at the highest subsonic Mach number. 

3. Optimize for the design parameters that minimize control surface activity, putting no active constraints on the 
values of ui n . 
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4. Introduce the minimum singular values into the target function (eq. (13)) and reoptimize the design variables, 
with initial values set to the optimal values obtained in step 2. 

5. Repeat step 4, if necessary, using different values for W a and ct d until either a desired value for a is obtained, 
or until a reaches an asymptotic value. 

6. Introduce the actuator and use as initial values (for the design parameters) the values associated with the 
optimized control law obtained in step 5, together with the actuator compensation terms, as shown earlier. 
Reoptimize to obtain the final optimum values for the design parameters. 

7. Test the performance of the resulting control law over a range of Mach number and flight configurations. 

This procedure not only yields the required control law, but it also provides the designer with physical insight 
regarding the system. Step 2, for example, can be skipped, although it provides the designer with the minimum 
value of control surface activity irrespective of robustness requirements. This is useful since, as already shown, 
robustness is achieved at the cost of increased control surface activity. Steps 3-6 can also be skipped. Step 3 can 
be replaced with a constrained optimization, with uj n limited between 0.65 up and 1.45 ujp if this method is to be 
used in a preliminary design stage in conjunction with a large structural-aeroservoelasticity optimization package. 
Steps 3-6 can be left to the final design stage in this latter case. 

Finally, in many instances the optimizer went to the absolute minimum of control surface activity associated 
with zero gains; that is, ai = ar = 0. This was easily avoided by imposing a lower bound on the gains. In this 
report, a? > 0 .5 yields satisfactory results, with the lower bound never turning to produce an active constraint. 

CONCLUDING REMARKS 

The aerodynamic energy approach yields low-order control laws with performances superior to those obtained 
using optimal control theory-based synthesis methods. Unlike previous works associated with the aerodynamic en- 
ergy method, there is no single specific control law, but instead an infinite number of control laws that result from 
the aerodynamic energy method. Also, higher order control laws, while still maintaining relatively low orders, pro- 
vide the system with more freedom to reduce control surface activities with a desired level of robustness. With 
sensors positioned around the locations required by the aerodynamic energy approach, all robust control laws, ir- 
respective of the method of synthesis, yield a large component of trailing-edge control rotation in phase with the 
velocity of its driving sensor. This is true over a wide range of frequencies that span the flutter frequency. The latter 
result is in accordance with the damping-type transfer function results obtained while developing the aerodynamic 
energy approach. 

A redefined synthesis procedure is suggested that replaces previous constraints imposed by the aerodynamic 
energy approach, with a target function that includes minimum singular values. The redefined procedure gives some 
weight to robustness while optimizing the design parameters. This is important since it is shown that robustness is 
achieved at a cost of increased control surface activity. 

It is hoped that additional examples relating to specific mathematical models will reconfirm the superiority of 
the aerodynamic energy approach regarding the synthesis of low-order control laws for flutter suppression. 
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APPENDIX A 

SUMMARY OF CONTROL LAWS USING THE INBOARD SENSOR 
Control Laws for System With No Actuator 

Control Law Al, optimized at Q D = 142.6 lb/ft 2 , w n > 2 , W, = 0 


S _ 4x2.1s 2 

SZ ~ s 2 + (2 x 2 x 0.90)s + (2) 2 


Control Law A2, optimized at Q D = 142.6 lb/ft 2 , w n > 35, W a = 0 

<5 _ 4 x 1 .2 s 2 

higb. ~ s 2 + (2 x 35 x 0.29)s + (35) 2 

Control Law A3, optimized at Q D = 142.6 Ib/ft 2 , w„ > 35 , = 1000 

6 4 x 1.03s 2 

5 2 + (2 x 35 x 0.34)s+ (35) 2 

Control Law A4, optimized at Q D = 124 )b/ft 2 , w„ > 35 , = 1000 


S 4 x 1.23s 2 

" s 2 + (2 X 39.3 x l)s+ (39 ,3) 2 

Control Law for System With Actuator and With Structural Filter 
Control Law AS, optimized at Q D = 142.6 Ib/ft 2 , W a = 1000 


Jc _ 4 x 1 x 630 (s + 91 .8) (s 2 + 270 ,9 s + 52,353) 

hgL 91.8 x 52,353 (s+ 630) [s 2 + (2 x 25.2 x 0.29)s+ (25 .2) 2 ] 

Control Law for System With Actuator But With No Structural Filter 
Control Law A6, optimized at Q D = 142.6 Ib/ft 2 , W a = 1000 


Sc _ 4 x 1.29 x 630 (s+ 186.6) (s 2 + 153.6s + 31,494) 

jW 186.6 x 31,494 (s+ 630) [s 2 + (2 x 15.8 x 0.39)s+ (15.8) 2 ] 


( A-l) 


(A-2) 


( A-3) 


( A-4) 


(A-5) 


(A-6) 
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APPENDIX B 

SUMMARY OF CONTROL LAWS USING THE OUTBOARD SENSOR 
Control Laws for System With No Actuator 


Control Law Bl, optimized at Qd = 142.6 lb/ft 2 , w n > 2 , W a - 0 


8 4 x 1 .72 s 2 

" s 2 + (2 x 2 x 0.89)s + (2) 2 


(B-l) 


Control Law B2, optimized at Qd = 142.6 lb/ft 2 , w n > 35 , W a - 0 


8 4 x 0 .96 s 2 

hngh = s 2 + (2 x 35 x 0.31)5+ (35) z 


(B-2) 


Control Law B3, optimized at Qd - 142.6 lb/ft 2 , w n > 35 , W 3 — 1000 


8 4 x 0 .69 s 2 

7 y = s 2 + (2 x 35 x 0.31)s+ (35) 2 


(B-3) 


Control Law B4, optimized at Qd = 124 lb/ft 2 , w n > 35 , W a — 1000 


8 4 x 0.72s 2 

7 y " s 2 + (2 x 35 x 0.98)s+ (35) 2 

Control Law for System With Actuator and With Structural Filter 
Control Law B5, optimized at Qd = 142.6 lb/ft 2 , W a = 1000 



4 x 0 .7 x 630 ( s + 126 .6) (s 2 + 251 ,8s + 42,602) 

126.6 x 42,602 ( 7 + 630) [s 2 + (2 x 28.7 x 0.26)s+ (28.7) 2 ] 


Control Laws for System With Actuator But With No Structural Filter 


Control Law B6, optimized at Qd = 142.6 lb/ft 2 , W a - 1000 


(B-4) 


(B-5) 



4 x 0.96 x 630 (s+ 185.8) (s 2 + 155s+ 29,658) 

185.8 x 29,658 (s+ 630) [s 2 + (2 x 16.5 x 0.41)s+ ( 16 .5) 2 ] 


(B-6) 
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Control Law B7, optimized at Q D = 151 lb/ft 2 , W, = 1000 


8 C _ 4 x 1.06 x 630 (s+ 189.3) (s 2 + 156 . 65 + 30,000) 

hputh ~ 189.3 x 30,000 (s+ 630) [s 2 + (2 x 15 x 0.37)s + ( 15) 2 ] 


(B-7) 
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Table 1. Comparisons between existing and reconstructed results for the DAST-ARW1 mathematical model 
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B3.Q = 124 lb/ft 2 3.03 166.2 - 8.88 14.2 -82 63 0.678 124 lb/ft 2 optimized at Q D = 142.6 lb/ft 2 . 

A4 2.66 147.5 -5.98 oo -58 80 0.834 Optimized at Q D = 124 lb/ft 2 . 

B4 2.67 150.2 -5.2 30.8 -63 67 0.82 ui. > 35. IV, = 1000. 
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Figure 2. Open-loop root locus plot at M = 0.9. 




Real part, rad/sec 


Figure 3. Root locus plot at M = 0.9; system off, arrows indicate increasing dynamic pressure. Reproduced from 
Newsom et al. (1980). 
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r 9018 

Figure 5. Qosed-loop root locus plot at M = 0.9; control law defined by equation (9), arrows indicate increasing 
dynamic pressure. Reproduced from Newsom ct al. (1980). 


MSV = 0.606 
A = -48 deg 



Imaginary, 

rad/sec 


Real, rad/sec 


9019 


Figure 6. Nyquist plot at Qd= 142.6 lb/ft 2 using control law given by equation (3), Newsom et al. 
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MSV = 0.64 

A = -6.0 dB gain margin 
B = 16.5 dB gain margin 
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’ 9020 

Figure 7. Nyquist plot of open-loop transfer function G( ico) H(iuj ) ; M = 0.9, Q = Q ma x ; H(iui) defined by 
equation (9), arrows indicate increasing frequency. Reproduced from Newsom et al. (1980). 
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Figure 9. Closed-loop dynamic pressure root locus plot of 65th-ordcr plant plus 4th-order control law. Designed 
with 25th-order plant, R u = 0.00001. Reproduced from Mukhopadhyay et al. (1981). 
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Figure 10. Nyquist plot at Qd = 151 lb/ft 2 using control law given by equation (30), Mukhopadhyay et al. (1981). 
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Figure 11. Nyquist plot of 65th-order plant model plus 4th-ordcr control law. Designed with 25th-order plant model, 
Q = 160 lb/flt 2 and R u = 0.00001. Reproduced from Mukhopadhyay et al. (1981). 
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Figure 1 2. Nyquist plot at Qo = 142.6 lb/ft 2 using aerodynamic energy control law A3. Optimized with > 35 
rad/sec, W, = 1000, and a D = 0.8 using inboard accelerometer. 
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Figure 13. Nyquist plot at Q D = 142.6 lb/ft 2 using aerodynamic energy control law B3. Optimized with > 35 
rad/sec, W s = 1000, and c r D = 0.8 using outboard accelerometer. 
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Figure 14. Root locus plot using aerodynamic energy control law A3. Optimized at Qd = 142.6 lb/ft 2 , with W s 
1000, g_ D = 0.8, and u > 35 using inboard accelerometer. 
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Figure 15. Root locus plot using aerodynamic energy control law B3. Optimized at Q& = 142.6 lb/ft 2 , with W s = 
1000, g_ D = 0.8, and > 35 using outboard accelerometer. 
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Figure 16. Nyquist plot at Qd = 142.6 lb/ft 2 using aerodynamic energy control law A5. Optimized with actuator 
and structural filter at Qd = 142.6 lb/ft 2 , with W a = 1000, g_ D = 0.9, and > 25 using inboard accelerometer. 
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Figure 17. Nyquist plot at Qd = 142.6 lb/ft 2 using aerodynamic energy control law B5. Optimized with actuator 
and with structural filter at Qd — 142.6 lb/ft 2 , with W s - 0, o_d = 0.85 using outboard accelerometer. 
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Figure 18. Nyquist plot at Q = 134 lb/ft 2 of aerodynamic energy control law A5. Optimized at Qo 
with actuator and structural filter using inboard accelerometer. 
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MSV = 0.682 
A = - 119 deg 
B = 63 deg 
C = - 9.04 dB 
D = 9.92 dB 
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Figure 19. Nyquist plot at Q = 124 lb/ft 2 of ae 
with actuator and structural filter using inboard 
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energy control law A5. Optimized at Qd = 142.6 lb/ft' 
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Figure 20. Nyquist plot at Q = 134 lb/ft 2 of aerodynamic energy control law B5. Optimized at Qd 
142.6 lb/ft 2 , with actuator and structural filter using outboard accelerometer. 
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Figure 21. Nyquist plot at Q = 124 lb/ft 2 of aerodynamic energy control law B5. Optimized at Qd 
142.6 lb/ft 2 , with actuator and structural filter using outboard accelerometer. 



Figure 22. Closed-loop root locus plot using aerodynamic energy control law A5. 
142.6 lb/ft 2 , with actuator and with structural filter at W a = 1000, o D = 0.9, and 
accelerometer. 
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Figure 23. Root locus plot using aerodynamic eneigy control law B5. Optimized at Q D = 142.6 lb/ft 2 , with actuator 
and structural filter at W 3 = 1000, g_ D = 0.85, and w-r > 25 using outboard accelerometer. 
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Figure 24. Nyquist plot at Qu = 151 lb/ft 2 using aerodynamic energy control law B7. Optimized with actuator but 
without structural filter at Q o = 151 lb/ft 2 , with W, = 1000 and g_ D = 1 using outboard accelerometer. 
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Figure 25. Closed-loop root locus plot using aerodynamic energy control law B8. Optimized at Q D = 151 lb/ft' 
with actuator but without structural filter using outboard accelerometer. 
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(a) Imaginary part. 

Figure 26. Real and imaginary parts of control deflection per unit displacement of outboard sensor using control 
law given by equation (9), Newsom et al. (1980). 
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Figure 26 Concluded. 
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(a) Imaginary part. 



Frequency, rad/sec ^ 


(b) Real part. 

Figure 27. Real and imaginary parts of control deflection per unit displacement of outboard sensor using control 
law given by equation (30), Mukhopadhyay et al. (1981). 
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(a) Imaginary part. 
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(b) Real part. 

Figure 28. Real and imaginary parts of control deflection per unit displacement of outboard sensor using control 
law given by equation (29), Mukhopadhyay et al. (1981). 
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(b) Real part. 

Figure 29. Real and imaginary parts of control deflection per unit displacement of inboard sensor using control 
law A3 with no actuator. 
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(a) Imaginary part. 
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(b) Real part. 

Figure 30. Real and imaginary parts of control deflection per unit displacement of outboard sensor using control 
law B3 with no actuator. 
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(b) Real part. 

Figure 31. Real and imaginary parts of control deflection per unit displacement of inboard sensor using control 
law A5. Optimized with actuator and with structural filter. 
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(a) Imaginary part. 
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(b) Real part. 

Figure 32. Real and imaginary parts of control deflection per unit displacement of outboard sensor using control 
law B5. Optimized with actuator and with structural filter. 
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(a) Imaginary part. 
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(b) Real part. 

Figure 33. Real and imaginary parts of control deflection per unit displacement of outboard sensor using control 
law B7. Optimized at Q D = 15 1 lb/ft 2 , with actuator but without structural filter. 
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